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Global  Modeling 

of  Internal  Tides 


Within  an  Eddying  Ocean  General  Circulation  Model 


BY  BRIAN  K.  ARBIC, 
PATRICK  C.  TIMKO,  E.  JOSEPH 


JAMES  G.  RICHMAN,  JAY  F.  SHRIVER, 
METZGER,  AND  ALAN  J.  WALLCRAFT 


ABSTRACT.  Ocean  tides,  and  the  atmospherically  forced  oceanic  general 
circulation  and  its  associated  mesoscale  eddy  field,  have  long  been  run  separately 
in  high-resolution  global  models.  They  are  now  being  simulated  concurrently  in  a 
high -resolution  version  of  the  HYbrid  Coordinate  Ocean  Model  (HYCOM).  The 
incorporation  of  horizontally  varying  stratification  with  the  addition  of  atmospheric 
forcing  yields  internal  tides  (internal  waves  of  tidal  frequency)  in  high-latitude, 
low- stratification  regions  that  are  qualitatively  different  from  those  in  earlier  global 
internal  tide  models,  in  which  atmospheric  forcing  and  horizontally  variable 
stratification  were  absent.  The  internal  tides  in  the  new  concurrent  HYCOM 
simulations  compare  well  with  those  measured  in  along-track  satellite  altimeter  data. 
The  new  concurrent  simulations  demonstrate  that  the  wavenumber  spectrum  of 
sea  surface  height — a  measure  of  the  energy  contained  in  different  length  scales — is 
dominated  in  some  locations  by  internal  tides  and  in  others  by  mesoscale  eddies. 
Tidal  kinetic  energies  in  the  new  concurrent  simulations  compare  well  with  those 
in  current-meter  observations,  as  long  as  sufficient  spatial  averaging  is  performed. 
The  new  concurrent  simulations  are  being  used  in  the  planning  of  future-generation 
satellite  altimeters,  in  the  provision  of  boundary  conditions  for  coastal  ocean  models, 
and  in  studies  of  ocean  mixing. 
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TIDES  AND  GLOBAL 
TIDE  MODELING 
Ocean  tides  result  from  the  difference 
in  the  gravitational  potential  of  the 
Moon  and  Sun  across  Earth’s  extent 
(Cartwright,  1999).  The  tidal  potential, 
or  equilibrium  tide,  has  a  simple  shape 
first  elucidated  by  Newton.  Actual  ocean 
tides  represent  a  complex  dynamical 
response,  first  elucidated  in  mathemati¬ 
cal  form  by  Laplace,  to  the  equilibrium 
tidal  forcing.  Ocean  tides  are  affected  by 
several  factors,  including  Earth’s  rota¬ 
tion,  friction  with  the  seafloor,  water 
depth  (which  controls  the  speed  of 
shallow- water  gravity  waves),  and  the 
presence  of  continents,  as  well  as  seem¬ 
ingly  exotic  effects  such  as  tides  in  the 
solid  Earth  and  perturbations  in  the 
gravitational  self-attraction  of  mass  in 
the  ocean  and  solid  Earth  brought  about 
by  tidal  motions  (Hendershott,  1972). 

The  shallow- water  equations  used  by 
modern  scientists  to  model  ocean  tides 
are  essentially  updated  versions  of  the 
Laplace  tidal  equations,  in  which  the 
effects  described  above  are  included.  An 
important  assumption  in  the  shallow- 
water  equations  is  that  motions  have 
horizontal  scales  that  are  large  compared 
to  the  water  depth  and  therefore  lie  in 
hydrostatic  balance  (Pedlosky,  1987; 
Vallis,  2006). 

Early  global  models  of  the  tides 
(e.g.,  Hendershott,  1972)  were  baro- 
tropic,  meaning  that  the  relatively  small 
variations  in  ocean  density  from  the  top 
to  the  bottom  of  the  water  column  were 
ignored.  In  a  shallow-water  model,  the 
assumption  of  constant  density  implies 
that  the  horizontal  velocities  do  not 
change  with  vertical  position  in  the 
water  column,  though  they  do  change 
with  latitude  and  longitude.  The  spacing 


between  adjacent  horizontal  grid  points 
in  early  global  tide  models  was  about 
6°  (670  km).  This  relatively  coarse  grid 
spacing  was  dictated  by  the  limited 
computer  power  available  at  the  time. 
Modern  computer  power  allows  compu¬ 
tation  of  the  global  tides  for  grid  spacing 
as  fine  as  l/25°-l/12°  (4.4-8. 9  km). 

As  has  been  recognized  for  some 
time  (e.g.,  Wunsch,  1975),  tidal  flow 
over  topographic  features  generates  an 
internal  tide.  In  a  stratified  fluid,  the 
vertical  motions  induced  by  flow  over 
topography  take  the  form  of  internal 
waves  (waves  along  interfaces  between 
fluids  of  differing  densities),  which 
are  the  subject  of  this  special  issue  of 
Oceanography.  Internal  tides  are  internal 
waves  with  tidal  frequencies.  Internal 
tides  can  have  displacement  amplitudes 
greater  than  50  m  and  current  speeds 
greater  than  2  m  s_1.  Internal  tides  are 
a  topic  of  great  current  interest,  for  a 
number  of  reasons.  The  vertical  motions 
(displacements)  associated  with  internal 
tides  affect  the  motions  of  gliders,  floats, 
and  other  oceanographic  instruments. 
They  also  impact  the  motions  of  subma¬ 
rines  and  acoustic  (sound)  waves  in  the 
ocean.  Most  importantly,  the  breaking 
of  internal  waves  is  a  leading  source  of 
mixing,  and  internal  tides  are  thought 
to  be  one  of  the  largest  energy  sources 


powering  mixing  in  the  deep  ocean. 

It  is  therefore  thought  that  tides  exert 
important  controls  on  the  stratification 
and  circulation  of  the  deep  ocean  (Munk 
and  Wunsch,  1998). 

Internal  tides  have  much  smaller  hori¬ 
zontal  scales  than  barotropic  tides,  and 
by  definition  can  only  be  numerically 
modeled  in  a  system  with  more  than  one 
layer  in  the  vertical  direction.  Because 
both  horizontal  and  vertical  resolution 
must  be  increased,  the  computational 
effort  required  to  model  internal  tides  is 
commensurately  greater  than  the  effort 
required  to  model  barotropic  tides.  One 
approach  to  this  challenge  is  to  model 
internal  tides  on  a  regional  scale,  which 
permits  higher  spatial  resolution  than  is 
possible  in  global  models.  Carter  et  al. 
(2012,  in  this  issue)  provide  a  review  of 
regional  internal  tide  models.  Because 
of  the  great  computational  expense 
involved  in  global  internal  tide  models, 
the  first  papers  on  global  modeling  of 
internal  tides  are  less  than  a  decade  old, 
and  the  number  of  published  papers 
that  focus  on,  or  at  least  include,  global 
models  of  internal  tides  is,  to  the  best  of 
our  knowledge,  still  quite  small  (Arbic 
et  al.,  2004;  Simmons  et  al.,  2004;  Hibiya 
et  al.,  2006;  Simmons,  2008;  Arbic  et  al., 
2010;  Niwa  and  Hibiya,  2011).  In  the 
earliest  global  internal  tide  modeling 
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studies,  astronomical  tidal  forcing  was 
the  only  forcing  present  in  the  model. 
Due  to  the  lack  of  atmospheric  forcing, 
which  sets  the  oceanic  stratification,  the 
oceanic  stratification  was  assumed  to  be 
horizontally  uniform.  Oceanic  stratifica¬ 
tion  in  subtropical  regions  was  assumed 
to  hold  everywhere  in  the  ocean  and, 
as  a  result,  the  internal  tides  in  polar 
regions  in  early  global  internal  tide 
models  were  noted  to  be  quite  inaccurate 
(Padman  et  al.,  2006). 

As  in  many  subdisciplines  of  physi¬ 
cal  oceanography,  the  study  of  tides  has 
been  revolutionized  by  the  advent  of 
satellite  altimetry  (Fu  and  Cazenave, 
2001).  From  an  orbit  of  approximately 


1,300  km,  satellite  altimeters  can  mea¬ 
sure  sea  surface  heights  averaged  over  an 
area  of  several  square  kilometers  with  an 
accuracy  of  about  1  cm.  This  accuracy 
implies  that  any  oceanographic  phenom¬ 
enon  with  a  sea  surface  height  (SSH)  sig¬ 
nature,  including  tides,  can  be  examined 
in  satellite  altimeter  data.  Global  baro- 
tropic  tide  models  that  are  developed 
from  satellite  altimeter  data  (e.g.,  Egbert 
et  al.,  1994;  Ray,  1999)  compare 
extremely  well  to  independent  observa¬ 
tions  such  as  tide  gauges  and  bottom 
pressure  recorders  (Shum  et  al.,  1997). 
Because  tides  contribute  about  80%  of 
the  SSH  variance  measured  by  altim¬ 
eters,  examination  of  nontidal  oceanic 
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Figure  1.  (a)  Amplitude  (cm)  of  M2  surface  tidal  elevation  in  TPXO  (Egbert  et  at,  1994),  a 
highly  accurate  barotropic  tide  model  constrained  by  satellite  altimetry,  (b)  Amplitude  (cm) 
of  M2  surface  tidal  elevation  in  HYbrid  Coordinate  Ocean  Model  (HYCOM)  simulations  in 
which  the  tide  is  unconstrained  by  satellite  altimetry.  White  lines  in  (a)  and  (b)  indicate  lines 
of  constant  phase.  Insets  to  the  right  of  (a)  and  (b)  display  amplitude  and  phase  in  the  west¬ 
ern  Pacific  region  delineated  by  boxes  in  (a)  and  (b).  (c)  January  19,  2012,  snapshot  of  sea 
surface  height  (cm)  in  the  real-time  data-assimilative  HYCOM  without  tides. 


motions  in  altimeter  data  can  only 
take  place  after  tides  have  been  accu¬ 
rately  removed  from  altimeter  records 
(Le  Provost,  2001).  This  illustrates  a  fun¬ 
damental  truth  about  the  place  of  tides 
in  physical  oceanography:  they  are  “sig¬ 
nal”  in  many  applications,  and  “noise”  in 
many  other  applications.  Satellite  altim¬ 
eters  have  also  been  instrumental  in  the 
study  of  internal  tides.  Although  internal 
tides  take  on  their  greatest  amplitudes  at 
depth,  their  signature  in  altimeter  data  is 
large  enough  to  be  detectable  (e.g.,  Ray 
and  Mitchum,  1996, 1997). 

Barotropic  tides  generate  internal 
tides,  and  internal  tides  in  turn  feed  back 
onto  the  barotropic  tides.  Inferences 
from  altimetry-constrained  barotropic 
tide  models  show  that  about  one-third 
of  global  tidal  energy  dissipation  occurs 
in  regions  of  rough  topography,  where 
internal  tides  are  generated  (Egbert  and 
Ray,  2000).  Internal  tide  generation  thus 
acts  as  a  damping  mechanism  for  the 
barotropic  tides. 

Figure  la,b  shows  the  surface  eleva¬ 
tion  amplitude  of  the  principal  lunar 
semidiurnal  tide  M2,  the  largest  tidal 
constituent  in  the  ocean.  Figure  la  dis¬ 
plays  the  M2  amplitude  in  TPXO  (Egbert 
et  al.,  1994),  a  highly  accurate  altimetry- 
constrained  barotropic  tide  model. 

Figure  lb  displays  the  M2  amplitude  in 
the  new  concurrent  HYbrid  Coordinate 
Ocean  Model  (HYCOM)  simulations 
described  in  this  paper.  Although  the 
HYCOM  M2  amplitudes  are  clearly  very 
similar  to  those  in  TPXO,  they  are  less 
accurate  due  to  the  lack  of  constraints 
applied  using  satellite  altimetry.  (Thus 
far,  our  HYCOM  simulations  are  “for¬ 
ward”  simulations,  that  is,  simulations 
unconstrained  by  data.)  Another  key 
difference  is  that  the  HYCOM  simula¬ 
tions  include  internal  tides,  which  can  be 
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seen  in  Figure  lb  as  small-scale  “wiggles” 
in  the  lines  of  constant  phase  (white 
lines  in  Figure  la,b)  and  as  small-scale 
“dimples”  in  the  amplitude  contours  of 
Figure  lb.  The  presence  of  internal  tide 
perturbations  to  amplitude  and  phase  in 
the  new  concurrent  HYCOM  simulation, 
and  the  lack  of  these  perturbations  in 
TPXO,  is  more  easily  seen  by  contrast¬ 
ing  the  insets  to  Figure  la,b.  The  insets 
focus  on  the  western  Pacific,  a  region  of 
strong  internal  tides. 

GLOBAL  GENERAL 
CIRCULATION  MODELING 
With  the  notable  exception  of  gravita¬ 
tionally  forced  ocean  tides,  the  atmo¬ 
sphere  provides  the  primary  forcing 
of  large-scale  motions  in  the  ocean 
(motions  on  scales  from  1  km  up  to 
basin  scales).  For  instance,  exchanges 
between  the  ocean  and  atmosphere  at 
the  sea  surface  exert  a  strong  control  on 
ocean  stratification.  Seawater  is  strati¬ 
fied  by  density,  and  seawater  density  is 
a  function  of  temperature,  salinity,  and 
pressure.  Temperatures  in  the  upper 
ocean  (and  lower  atmosphere)  are  deter¬ 
mined  by  fluxes  of  sensible  and  latent 
heat  between  the  atmosphere  and  ocean 
(Yu  and  Weller,  2007).  Salinities  in  the 
upper  ocean  are  set  by  the  difference 
between  evaporation  and  precipitation 
at  the  ocean  surface  (Yu,  2007;  Schmitt, 
2008).  Because  the  buoyancy  (density) 
of  seawater  at  the  ocean  surface  is  con¬ 
trolled  by  temperature  and  salinity, 
the  forcing  described  above  is  known 
as  buoyancy  forcing.  The  atmosphere 
forces  the  ocean  through  pressure 
loading  (Ponte,  1993)  and  wind  stress 
(Pond  and  Pickard,  1983)  in  addition 
to  the  buoyancy  forcing.  Wind  stress  at 
the  ocean  surface  drives  near-inertial 
motions  (Alford,  2003;  Simmons  and 


Alford,  2012,  in  this  issue),  Ekman  flows 
(Pedlosky,  1996),  and  the  large-scale 
general  circulation  (Pedlosky,  1996). 
Winds  blowing  on  the  sea  surface  also 
drive  oceanic  motions  on  small  scales 
(scales  from  about  1  cm  to  about  1  km), 
such  as  capillary  waves  (Cox,  1958), 
surface  wind  waves  (Phillips,  1966),  and 
Langmuir  cells  (Langmuir,  1938).  These 
small-scale  motions,  however,  are  not 
resolved  in  models  such  as  HYCOM, 
which  make  the  hydrostatic  approxima¬ 
tion  (Pedlosky,  1987;  Vallis,  2006). 

The  oceanic  general  circulation  con¬ 
sists  of  the  well-known  current  systems 
around  the  globe,  such  as  the  Gulf 
Stream  in  the  western  North  Atlantic 
Ocean,  the  Kuroshio  in  the  western 
North  Pacific  Ocean,  the  Agulhas 
Current  off  of  Southern  Africa,  and  the 
Antarctic  Circumpolar  Current  in  the 
Southern  Ocean  (Schmitz,  1996a,b; 
Siedler  et  al.,  2001).  On  timescales  of 
about  10-200  days,  these  currents  mean¬ 
der  and  generate  highly  energetic  meso- 
scale  eddies  (Schmitz,  1996a,b;  Stammer, 
1997),  the  spinning  oceanic  dynamical 
counterparts  of  atmospheric  weather 
systems  generated  through  a  hydro- 
dynamical  process  known  as  baroclinic 
instability  (Pedlosky,  1987;  Vallis,  2006). 
Because  currents  meander  and  meso- 
scale  eddies  are  born  and  die  over  rela¬ 
tively  long  timescales  (scales  much  lon¬ 
ger  than,  for  instance,  tidal  timescales), 
these  currents  and  eddies  are  sometimes 
referred  to  as  low-frequency  motions. 

A  substantial  percentage  of  the  kinetic 
energy  in  the  ocean  resides  in  mesoscale 
eddies  and  other  transient  low-frequency 
features  such  as  current  meanders. 
Because  the  horizontal  length  scales  of 
strong  currents  and  mesoscale  eddies  are 
relatively  small  (about  100  km),  numeri¬ 
cal  models  of  the  ocean  must  have  high 


resolution  in  order  for  the  modeled 
low-frequency  kinetic  energies  to  com¬ 
pare  well  to  those  recorded  in  observa¬ 
tions  from  satellite  altimeters  and  other 
instruments  (Hurlburt  and  Hogan,  2000; 
Maltrud  and  McClean,  2005).  The  grid 
spacing  in  current  state-of-the-art  global 
high-resolution  ocean  models  is  about 
1/10°  to  1/12°,  which  is  thought  to  be 
sufficient  to  capture  most  of  the  oceanic 
mesoscale.  Indeed,  ocean  models  having 
resolutions  this  high  or  higher  are  known 
as  “eddying”  models  (Hecht  and  Hasumi, 
2008).  Note,  however,  that  Ihoppil  et  al. 
(2011)  find  that  increasing  the  resolution 
of  HYCOM  from  1/12°  to  1/25°  yields 
further  increases  in  the  energy  of  mod¬ 
eled  low-frequency  flows,  suggesting 
that  resolutions  of  1/25°  or  higher  might 
be  necessary  to  fully  resolve  mesoscale 
eddies  and  current  meanders. 

HYCOM  is  currently  being  devel¬ 
oped  for  use  as  the  next-generation 
global  ocean  forecast  model  for  the 
United  States  Navy  (Chassignet  et  al., 
2007;  Metzger  et  al.,  2010;  http://hycom. 
org;  http://www7320.nrlssc.navy.mil/ 
GLBhycoml-12).  The  Navy  needs  to 
know  the  state  of  the  ocean — the  strength 
of  currents,  sea  surface  height,  clarity 
of  the  water,  and  sound  propagation 
properties,  among  others — on  a  range 
of  temporal  and  spatial  scales  at  diverse 
locations  around  the  globe.  The  model 
is  forced  by  fluxes  obtained  from  the 
Navy  Operational  Global  Atmospheric 
Prediction  System  (NOGAPS;  Rosmond 
et  al.,  2002).  The  data-assimilative  ver¬ 
sion  of  the  model  uses  satellite  altimet- 
ric  sea  surface  height,  radiometric  sea 
surface  temperature,  and  in  situ  ship, 
float,  and  profile  temperature  and  salin¬ 
ity  data  using  the  Navy  Coupled  Ocean 
Data  Assimilation  system  (NCODA; 
Cummings,  2005)  to  provide  a  statistical 
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blending  of  model  and  observations.  The 
global  ocean  model  runs  on  a  combina¬ 
tion  Mercator  and  combination  tripolar 
grid  (Murray,  1996)  and  Mercator  grid 
with  an  equatorial  grid  resolution  of 
0.08°  (8.9  km)  and  is  coupled  to  a  sea  ice 
model  at  high  latitudes.  HYCOM  simu¬ 
lations  at  0.04°  (4.5  km)  have  already 
begun  to  be  developed  and  run.  The 
model  reproduces  the  general  circula¬ 
tion  of  the  global  ocean,  the  strength 
and  variability  of  western  boundary 
currents  such  as  the  Gulf  Stream  and 
the  Kuroshio,  and  the  mesoscale  eddies 
generated  by  instabilities  of  the  major 
currents.  The  model  is  validated  by  com¬ 
parison  to  in  situ  observations  not  used 
in  assimilation  (Hurlburt  et  al.,  2010; 
Ihoppil  et  al.,  2011).  The  global  model 
performance  is  not  as  good  in  coastal 
regions  where  smaller-scale  features  in 
bathymetry,  and  new  phenomena  such  as 
tides,  become  important  in  determining 
the  state  of  the  ocean. 

To  contrast  nontidal  and  tidal  SSH 
fields,  Figure  lc  displays  a  snapshot 
of  SSH  in  the  data  assimilative  (ocean 
forecasting)  version  of  HYCOM,  which 
currently  does  not  include  tides.  The 
relatively  high  sea  surface  heights  in  the 
subtropical  regions  are  the  subtropical 
gyres  (Pedlosky,  1996).  The  strongest 
gradients  in  sea  surface  height  take 
place  in  western  boundary  currents  such 
as  the  Gulf  Stream  and  the  Kuroshio. 

The  numerous  small-scale  features 
present  in  Figure  lc  are  the  oceanic 
mesoscale  eddies. 

MARRIAGE  OF  GLOBAL  TIDAL 
AND  GENERAL  CIRCULATION 
MODELING 

Recently,  we  have  shown  that  tides 
and  the  atmospherically  forced  eddy¬ 
ing  ocean  general  circulation  can  be 


modeled  concurrently  on  a  global  scale; 
hence,  tides  and  nontidal  motions  now 
interact  in  the  model,  as  in  the  ocean. 
We  add  astronomical  tidal  forcing  of 
the  four  largest  semidiurnal  constitu¬ 
ents  (M2,  S2,  N2,  and  K2)  and  the  four 
largest  diurnal  constituents  (K1;  Op  P , , 
and  Qj)  to  the  HYCOM  general  circu¬ 
lation  model,  which  is  also  forced  by 
atmospheric  fields  as  discussed  in  the 
previous  section.  Arbic  et  al.  (2010) 
offer  more  details  of  the  new  concur¬ 
rent  tidally  and  atmospherically  forced 
simulations.  Here,  we  note  that  we  use 
the  parameterized  topographic  internal 
wave  drag  scheme  of  Garner  (2005), 
modified  as  described  in  Arbic  et  al. 
(2010),  in  order  to  limit  the  maximum 
decay  rates  to  (9  hrs)-1.  The  wave  drag 
scheme  parameterizes  the  drag  on  tidal 
flows  resulting  from  the  generation  of 
unresolved  small-vertical-scale  internal 
waves  by  tidal  flow  over  rough  topog¬ 
raphy.  Because  the  wave  drag  acting  on 
tidal  motions  is  known  to  have  a  differ¬ 
ent  strength  than  the  wave  drag  acting 
on  nontidal  motions  (Bell,  1975),  we 
separate  the  tidal  and  nontidal  bottom 
flows  in  our  model  using  a  25-hour  run¬ 
ning  boxcar  filter.  The  self-attraction 
and  loading  term  (Hendershott,  1972)  is 
simplified  using  the  scalar  approxima¬ 
tion  (Ray,  1998)  and  is  applied  only  to 
the  nonsteric  component  of  sea  surface 
height,  which  is  dominated  by  the  large- 
scale  barotropic  tides. 

With  these  adaptations,  HYCOM 
generates  both  barotropic  and  internal 
tides  amidst  the  eddying  general  circu¬ 
lation.  The  barotropic  tides  in  the  new 
concurrent  HYCOM  tidal  simulations 
were  compared  against  the  standard  set 
of  102  pelagic  tide  gauges  (Shum  et  al., 
1997)  in  Arbic  et  al.  (2010),  who  found 
that  HYCOM  captured  about  93%  of 


the  sea  surface  elevation  variance  of  the 
eight  largest  tidal  constituents  in  the 
tide  gauge  records.  Arbic  et  al.  (2010; 
see  their  Figure  3)  showed  that  the 
nontidal  sea  surface  height  variability 
was  not  adversely  affected  by  the  pres¬ 
ence  of  tides  in  the  model.  In  the  new 
concurrent  HYCOM  tidal  simulations, 
in  contrast  to  early  internal  tide  simula¬ 
tions  (Arbic  et  al.,  2004;  Simmons  et  al., 
2004),  the  atmospheric  forcing  present 
permits  a  realistic  horizontally  varying 
stratification  in  the  ocean  model.  As  a 
result,  the  internal  tides  in  high  latitudes 
are  generated  in  a  stratification  that  is 
more  realistic  than  was  possible  in  early 
internal  tide  simulations.  In  Figure  2,  we 
display  the  consequences  of  this  horizon¬ 
tally  varying  stratification  for  M2.  The 
figure  displays  the  amplitude  of  the  M2 
internal  tide  signature  in  the  steric  sea 
surface  height  (i.e.,  the  sea  surface  height 
component  having  to  do  with  stratifica¬ 
tion  effects).  Figure  2a,  taken  from  a 
HYCOM  simulation  like  those  in  earlier 
global  internal  tide  studies — using  a  two- 
layer,  horizontally  uniform  stratification 
typical  of  the  subtropics — displays  an 
artificially  strong  internal  tide  in  high- 
latitude  regions  such  as  the  Labrador  Sea 
and  Drake  Passage.  In  the  new  concur¬ 
rent  HYCOM  simulations  with  tides 
embedded  within  an  atmospherically 
forced  model  (Figure  2b),  this  defect  is 
removed.  Figure  2b  also  demonstrates 
that  the  introduction  of  a  more  realistic 
32-layer  horizontally  varying  stratifica¬ 
tion  leads  to  increased  internal  tide 
amplitudes  at  most  of  the  major  genera¬ 
tion  sites  in  the  mid-  and  low-latitude 
oceans,  when  compared  to  the  internal 
tides  generated  in  a  two-layer  hori¬ 
zontally  uniform  stratification.  See,  for 
instance,  the  area  around  Hawai’ i. 

Using  the  new  concurrent  tide  model, 
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we  have  begun  to  examine  a  host  of  sci¬ 
entific  and  operational  questions  about 
internal  tides.  In  the  following  sections, 
we  give  brief  overviews  of  some  of  our 
preliminary  results  of  this  exploration. 


averaged  the  root-mean-square  (rms) 
perturbation  magnitudes  over  the  five 
hotspot  regions  delineated  by  boxes  in 
Figure  3b,  and  found  that  the  modeled 
rms  values  agree  with  the  altimeter  rms 


values  to  within  about  20%  over  these 
regions.  All  of  this  supports  further 
confidence  in  the  realism  of  the  inter¬ 
nal  tides  in  our  HYCOM  simulations, 
in  the  same  manner  that  side-by-side 


GLOBAL  COMPARISON  OF 
MODELED  INTERNAL  TIDES 
WITH  SATELLITE  ALTIMETER 
ESTIMATES 

We  have  begun  to  validate  the  internal 
tides  in  our  model  with  a  global  compar¬ 
ison  to  satellite  altimeter  observations. 

As  in  the  previous  section,  we  focus  here 
on  M2,  the  largest  constituent  of  oceanic 
tides.  We  apply  a  50-400  km  band-pass 
filter  to  the  M2  amplitudes  at  the  sea 
surface  in  order  to  remove  both  small- 
horizontal-scale  noise  (in  the  altimeter 
data)  and  large-scale  barotropic  tides  (in 
the  altimeter  data  and  in  the  HYCOM 
simulation).  The  band  passing  there¬ 
fore  leaves  only  the  low-mode  internal 
tide  remaining.  Figure  3  displays  the 
results  of  applying  the  band-pass  filter 
to  the  surface  M2  elevations.  Figure  3a 
displays  results  from  along-track  satel¬ 
lite  altimeter  data  (Ray  and  Mitchum, 
1996,  1997;  Richard  Ray,  NASA,  pers. 
comm.,  2011),  while  Figure  3b  displays 
results  from  the  model  output  interpo¬ 
lated  to  the  same  altimeter  tracks.  The 
“hotspots”  of  internal  tide  generation  at 
prominent  topographic  features  such  as 
Hawai’ i  and  the  Aleutians  in  the  North 
Pacific,  the  Tuamotu  Archipelago  in  the 
South  Pacific,  Madagascar  in  the  Indian 
Ocean,  and  so  on,  are  very  clearly  seen 
in  both  panels  of  Figure  3.  The  internal 
waves  radiate  away  from  the  hotspots 
as  focused  beams,  which  propagate  for 
thousands  of  kilometers.  The  magnitudes 
of  the  internal  tide  perturbations  to  M2 
sea  surface  elevations  are  typically  about 
1  cm  or  more  in  the  hotspot  regions.  We 
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Figure  2.  Amplitude 
(cm)  of  the  M2  internal 
tide  signature  in  steric 
sea  surface  height  of 
HYCOM  experiments 
(a)  run  as  in  early  inter¬ 
nal  tide  simulations 
with  M2  tidal  forcing  as 
the  only  forcing  present 
and  with  a  horizontally 
uniform  two-layer 
stratification,  (b)  run 
with  astronomical  tidal 
forcing  present  amongst 
an  atmospherically 
forced  32-layer  general 
circulation  model  with 
horizontally  nonuniform 
stratification.  Boxes  in 
Figure  2b  denote  regions 
over  which  wavenumber 
spectra  are  computed 
in  Figure  4. 
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Figure  3.  M2  internal 
tide  perturbation  (cm) 
to  sea  surface  elevation. 
Perturbations  are  deter¬ 
mined  via  band  passing 
the  total  M2  sea  surface 
elevation  amplitude. 

(a)  Results  from  analysis 
of  along-track  satellite 
altimeter  data  (Ray  and 
Mitchum  1996, 1997; 
Richard  Ray,  NASA,  pers. 
comm.,  2011).  (b)  Results 
from  high-resolution  sim¬ 
ulation  of  HYCOM.  The 
boxes  in  (b)  delineate 
“hotspot”  regions  over 
which  root-mean-square 
values  of  the  internal  tide 
perturbations  are  com¬ 
puted  from  both  altim¬ 
eter  data  and  HYCOM 
output  (see  text). 
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plots  of  the  low-frequency  sea  surface 
height  variance  in  high-resolution  mod¬ 
els  versus  satellite  altimeter  data  (see, 
for  instance,  Figure  12  in  Maltrud  and 
McClean,  2005)  lend  confidence  in  the 
ability  of  high-resolution  models  to 
simulate  mesoscale  eddies. 

TIDAL  AND  NONTIDAL 
CONTRIBUTIONS  TO  THE  SEA 
SURFACE  HEIGHT  SPECTRUM 
We  have  begun  to  use  our  model  to  dis¬ 
tinguish  between  the  high-frequency 
tidal  and  low-frequency  nontidal  con¬ 
tributions  to  quantities  of  interest  such 
as  the  wavenumber  spectrum  of  sea 
surface  height.  Because  the  satellite 
altimeter  repeats  tracks  every  10  days, 
high-frequency  motions  such  as  the  tides 
are  aliased  into  longer  periods  (Parke 


et  al.,  1987).  In  satellite  altimeter  data, 
therefore,  it  is  sometimes  difficult  to 
distinguish  between  tidal  motions  and 
lower- frequency  motions.  Because  our 
model  output  is  written  at  hourly  inter¬ 
vals,  we  can  easily  separate  low-  from 
high-frequency  signals  in  the  model. 

In  Figure  4  we  show  the  wavenumber 
spectrum  of  SSH  from  the  model.  The 
wavenumber  spectrum  measures  the 
energy  content  of  a  signal  as  a  func¬ 
tion  of  wavenumber,  which  is  inversely 
proportional  to  the  wavelength.  The 
slope  of  the  wavenumber  spectrum 
is  of  great  theoretical  interest,  as  it  is 
used  to  infer  the  dominant  dynamics 
of  low-frequency  flows  (e.g.,  Le  Traon 
et  al.,  2008;  Xu  and  Fu,  2011).  In  regions 
where  low-frequency  motions  are  more 
energetic  than  tidal  motions,  such  as 


Figure  4.  Wavenumber  spectrum  of 
total  (black),  high-frequency  (red) 
and  low-frequency  (blue)  sea  surface 
height,  where  two  days  is  the  divid¬ 
ing  point  between  low  and  high 
frequencies,  for  (a)  Kuroshio  region, 
a  region  of  very  strong  low-frequency 
motions,  and  (b)  region  north  of 
Hawaii,  in  which  low-frequency 
motions  are  weaker  than  M2  internal 
tidal  motions.  The  two  regions  are 
delineated  by  the  westernmost  and 
easternmost  boxes,  respectively, 
in  the  North  Pacific  in  Figure  2b. 
Spectra  are  computed  from  a 
HYCOM  simulation  containing  both 
tides  and  low-frequency  atmospheri¬ 
cally  forced  motions.  Extra  green 
lines  are  drawn  in  at  the  best-fit 
slopes  of  -4.38  and  -3.87  over  low- 
frequency  motions  in  the  same 
wavenumber  band  discussed  in 
Xu  and  Fu  (201 1 ).  As  discussed  in 
Le  Traon  et  al.  (2008)  and  Xu  and  Fu 
(201 1 ),  -1 1/3  and  -5  slopes  imply 
surface  quasi-geostrophic  and  quasi- 
geostrophic  dynamics,  respectively. 


Wavenumber  (cpm) 


the  Kuroshio  (see  westernmost  box 
in  Figure  2b),  low-frequency  motions 
dominate  the  spectrum  of  total  SSH 
(Figure  4a).  However,  in  regions  where 
internal  tides  are  energetic,  for  instance 
near  Hawai’i  (see  easternmost  box  in 
Figure  2b),  high-frequency  tidal  motions 
dominate  the  high  end  of  the  wave- 
number  spectrum  (Figure  4b).  The  latter 
observation  implies  that  in  data  taken  by 
the  planned  high-resolution  wide-swath 
satellite  altimeter  (Fu  and  Ferrari,  2008), 
internal  tides  will  have  to  be  removed 
very  accurately  before  low-frequency 
oceanic  motions  can  be  studied.  Tides 
have  always  been  an  important  source  of 
noise  in  altimeter  measurements  of  SSH. 
However,  internal  tides  will  become  a 
more  important  source  of  noise  for  the 
wide-swath  altimeter  mission  than  for 
previous  altimeter  missions.  This  issue  is 
due  to  the  much  higher  spatial  resolution 
of  the  wide-swath  mission,  which  will 
collect  data  at  horizontal  scales  where  the 
internal  tides  contain  substantial  energy. 

GLOBAL  COMPARISON  OF 
MODELED  TIDAL  CURRENTS 
WITH  CURRENT  METER 
OBSERVATIONS 

We  have  begun  to  compare  the  tidal  cur¬ 
rents  in  our  model  to  those  in  current 
meter  observations.  The  set  of  current 
meter  observations  we  use  is  built  upon 
the  current  meter  database  in  Scott  et  al. 
(2010),  which  includes  observations 
from  1,528  locations,  denoted  by  blue 
dots  in  Figure  5a.  A  total  of  4,934  instru¬ 
ments  are  spread  out  over  many  depths 
at  the  1,528  locations. 

Because  the  ocean  is  stratified,  tidal 
currents  have  vertical  structure,  just  as 
tidal  displacements  do.  In  Figure  5b,  we 
compare  the  vertical  structure  of  the 
M2  tidal  kinetic  energy  in  our  HYCOM 
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simulations  against  that  seen  in  the  cur¬ 
rent  meters.  The  time-averaged  kinetic 
energy  is  calculated  over  the  M2  tidal 
period  at  each  instrument  location  for 
both  model  output  and  current  meter 
data,  and  then  averaged  over  all  locations 
having  instruments  within  a  specified 
depth  range.  As  can  be  seen  in  Figure  5b, 
after  this  spatial  averaging  is  performed, 
the  model  captures  the  vertical  structure 
of  the  semidiurnal  M2  tide  very  well — 
having  less  than  a  25%  relative  error  in 
the  top  2,000  m.  We  have  completed  a 
much  more  comprehensive  model-data 
comparison  for  tidal  currents  in  the 
North  Atlantic  (recent  work  of  author 
Timko  and  colleagues)  in  which  the 


observations  have  been  divided  on  a 
regional  basis.  The  study  of  the  North 
Atlantic  shows  that  model  performance 
can  vary  significantly  from  one  region 
to  another.  The  cause  of  regional  differ¬ 
ences  in  model  performance  is  a  subject 
of  future  investigation. 

Model  diurnal  kinetic  energies  appear 
substantially  weaker  than  in  observa¬ 
tions  (not  shown).  The  cause  appears  to 
be  that  tuning  the  topographic  internal 
wave  drag  for  M2  as  we  have  done  in 
Arbic  et  al.  (2010)  and  here  implies  that 
the  diurnal  tides  in  the  model  are  over¬ 
damped.  Theoretical  considerations 
(Bell,  1975)  imply  that  the  wave  drag  on 
diurnal  motions  is  less  strong  than  for 


Kinetic  Energy 


Figure  5.  (a)  Horizontal  loca¬ 
tions  of  current  meters  used 
to  validate  the  tidal  currents  in 
the  HYCOM  simulations.  Units 
of  x-  and  y-axes  are  degrees. 

(b)  Vertical  distributions  of 
the  M2  tidal  kinetic  energy, 
spatially  averaged  over  depth- 
binned  instruments,  for  current 
meter  observations  (“obs.")  and 
HYCOM  simulations  ("model"). 


semidiurnal  motions.  We  are  currently 
working  to  implement  a  frequency- 
dependent  wave  drag  into  the  HYCOM 
simulations.  It  is  not  simple  to  do  this, 
however,  because  the  model  runs  in  the 
time  domain  whereas  semidiurnal  and 
diurnal  tides  are  most  easily  separated  in 
the  frequency  domain. 

IMPLICATIONS  AND  FUTURE 
DEVELOPMENTS 
The  development  of  a  model  with  baro- 
tropic  and  internal  tides  embedded 
within  an  eddying  ocean  general  circula¬ 
tion  has  opened  up  a  host  of  scientific 
questions  and  applications,  some  of 
which  are  described  in  earlier  sections 
of  this  paper.  As  discussed,  the  model 
is  expected  to  be  useful  in  planning  for 
the  next-generation  wide-swath  satel¬ 
lite  altimeter  (Fu  and  Ferrari,  2008). 
Indeed,  the  model  already  has  been 
used  to  estimate  the  relative  importance 
of  tidal  versus  nontidal  contributions 
to  the  wavenumber  spectrum  of  sea 
surface  height  in  the  regions  where  a 
planned  airborne  test  campaign  will  be 
undertaken  (Ernesto  Rodriguez,  NASA, 
2011, pers.  comm.).  Three-dimensional 
maps  of  tidal  displacements  and  cur¬ 
rents  are  potentially  useful  for  removing 
tidal  “noise”  from  data  taken  by  various 
oceanographic  instruments.  In  addition, 
the  new  concurrent  tide  model  is  being 
used  in  the  Climate  Process  Team  project 
“Collaborative  Research:  Representing 
Internal-Wave  Driven  Mixing  in  Global 
Ocean  Models,”  which  focuses  on 
improving  estimates  of  internal-wave- 
mediated  mixing  in  the  ocean.  The  proj¬ 
ect  involves  several  institutions,  Jennifer 
MacKinnon  of  the  Scripps  Institution  of 
Oceanography  leads  the  project,  and  the 
National  Science  Foundation  funds  it. 

The  new  concurrent  tide  model  is  also 
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expected  to  be  useful  for  regional  and 
coastal  modeling.  Regional  and  coastal 
models  often  include  tidal  forcing,  but 
this  forcing  is  generally  barotropic,  such 
that  internal  tides  are  not  included  in 
the  open  boundary  conditions.  Kelly  and 
Nash  (2010)  showed  that  the  response 
of  shelf  tides  to  barotropic  plus  inter¬ 
nal  (i.e.,  fully  three-dimensional)  tidal 
forcing  can  be  quite  different  from  the 
response  to  forcing  that  includes  only 
barotropic  tides.  The  Australian  North 
West  Shelf  is  a  region  where  a  low 
vertical-mode  internal  tide  is  expected 
to  propagate  onto  the  shelf  from  the 
open  ocean  (Sam  Kelly,  University  of 
Western  Australia,  pers.  comm.,  2012). 
Current  meter  observations  (Holloway 
et  al.,  2001)  show  propagation  of  a 
mode- 1  internal  tide  onto  the  shelf.  The 
global  model  at  4  km  resolution  under¬ 
estimates  the  strength  of  the  internal 
tide  compared  to  the  observations,  but 
shows  onshore  propagation  from  deep 
water  onto  the  shelf  consistent  with  the 
observations.  However,  a  regional  model 
forced  by  barotropic  tidal  boundary 
conditions  shows  onshore  propagation 
of  the  internal  tide  only  inshore  of  the 
critical  slope,  and  offshore  propaga¬ 
tion  over  the  outer  continental  slope,  in 
disagreement  with  the  observations.  By 
including  tidal  forcing  in  a  fully  three- 
dimensional  global  ocean  circulation 
model,  we  will  provide  an  internal  tide 
capability  everywhere,  and  allow  nested 
models  to  include  internal  tides  at  their 
open  boundaries. 

A  number  of  improvements  are 
needed  for  the  new  concurrent  tide 
model.  The  self-attraction  and  loading 
term  (Hendershott,  1972)  is  imple¬ 
mented  in  a  rather  crude  way,  using  the 
so-called  scalar  approximation  (Ray, 
1998),  and  it  needs  to  be  improved  upon 
in  future  versions  of  the  model.  The 


frequency  dependence  of  topographic 
wave  drag  (Bell,  1975)  is  problematic 
because  the  model  is  run  in  the  time 
domain,  not  the  frequency  domain. 
Perhaps  the  greatest  technical  challenge 
remaining  is  how  to  perform  data  assim¬ 
ilation  on  the  tidal  as  well  as  nontidal 
motions  within  the  model.  Successful 
data  assimilation,  once  performed,  will 
make  the  model  much  more  useful  for 
operational  purposes  of  the  United  States 
Navy  and  other  users. 
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